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<^ ' Abstract 

^Q . We estimate the dependence of v^ to v e conversion on parameters 

#13 and bep for several experimental facilities studying neutrino os- 

£L2 ' cillations. We use the S-Matrix theory to estimate z7 e disappearance 

^vq ■ and compare estimates based on an older theory being used to extract 

#13 from the Double Chooz, Daya Bay, and RENO data, to assist in 
extracting an accurate value for #13 from these projects. We use val- 

k3i ■ ues of #13 within known limits, and estimate the dependence of v^ - 

5_i , v e CP violation (CPV) probability on 5cp in order to suggest new 

experiments to measure CPV for neutrinos moving in matter. 



1 Introduction 

In our present work we study v^ to v e neutrino conversion, v e disappearance, 
and CP violation (CPV) measurements using the S-Matrix method for neu- 
trino oscillations. The study of CP violation is essential for understanding 
weak interactions. Almost half a century ago CP violation in weak interac- 
tions was found in the decay of K\ into 7r + +7r~ [lj and 27r° [2], with branching 
ratios of the order of .001. The decay K\ — > it + v + v is almost entirely CP 
violating [3] but requires accurate determination of the CKM matrix jl] and 



accurate measurements. See Ref[5] for a review of this experiment and ref- 
erences. There have many other studies of CP asymmetries in weak decays: 
see Ref[6] for a recent study of B radiative decay with references to earlier 
work on CP violation in various weak decays. 

In recent years there have been a number of experimental studies of neu- 
trino oscillations using neutrino beams from accelerators and reactors, and 
important objectives of these experiments are the measurements of v^ to 
v e conversion and CPV. The first study in our present work is an estimate 
the v^ to v e conversion probability using parameters for the baseline and 
energy corresponding to MiniBooNEjT], JHF-Kamioka [8], MINOS [9], and 
CHOOZ [lOj-Double Chooz[llJ, which are on-going projects, although the 
CHOOZ project does not have a beam of v^ neutrinos. 

There have been many recent studies of CP and T symmetries via neu- 
trino oscillations for future facilities, e.g., see RefsJTSJ [T3], which also give 
references to earlier publications, and the ISS report [2] on future neutrino 
facilities. The two main parameters of interest in the present work are Sep, 
which is essentially unknown, and the angle 6*13. One possible future facility 
for studying CPV and the 5qp parameter is the LBNE Project, where neu- 
trino beams produced at Fermilab would have a baseline of L ~ 1200 km, 
being detected with deep underground detectors [151 [IB]- With the methods 
used in the present work, described below, predictions of CPV with the base- 
line and energies of the LBNE Project have recently been made for 5qp from 
90 to degrees [IT]. 

Although the angle # 13 is not well known, there are limits on its value. 
Recently the T2K collaboration [TS] published limits, and found a best-fit 
value of 6> 13 ~11 degrees, or sin^ = S13 ~ 0.19, which is one of the values 
we use in the present work on neutrino transition probability in Sec. 2. 

The angle #13 will be measured by the Daya Bay experiment |19[ 120] in 
China, the Double Chooz project J2TJ [22] in France, and RENO [23] in Korea, 
via v e disappearance. A very recent result from the Daya Bay project [23] 
concludes that s 13 ~ 0.15, and in our study (Sec. 3) of P e disappearance we 
use the the Daya Bay parameters to test the theory. Our results will also be 
useful for the Double Chooz and RENO projects. 

Finally, in Sec 4, using the expected range of values for # 13 , we estimate 
CPV for fx — e neutrino oscillation for the entire range of 5cp to help in the 
planning for future CPV experiments. 

A major complication for the determination of T, CP, and CPT violation 
is the interaction of neutrinos with matter as they travel along the baseline. 
These matter effects have been studied by a number of theorists. See, e.g., 
Refs[25j [261 ET| . One main objective of the present research is to estimate 
matter effects for the MiniBooNE, JHF-Kamioka, MINOS, and CHOOZ fa- 
cilities, to help find the values of 6*13 and Sep- 



For the basic interactions, which are CPT invariant for local theories, CP 
and T violation have the same magnitude. With matter effects T and CP are 
not directly related. Our present research is an extension of our recent work 
on T reversal violation [28J. In that study we used the formalism of Ref[29j 
for v e -H- Vp TRV, and that of Ref[30] for v e — > v^ conversion probability 
to calculate the effects of neutrinos moving through matter. In the present 
work we use the notation and formalism of Jacobson and Ohlsson[3Tj. who 
studied possible matter effects for CPT violation. 

CP violation in the a — b sector is given by the transition probability, 
denoted by V[y a — \ h>t,), for a neutrino of flavor a to convert to a neutrino of 
flavor b; and similarly for antineutrinos u a , v\>. The CPV probability differ- 
ences are defined as 

^ab P = V(v a ^V h )-V{v a ^V h ). (1) 

In our present work we study Viy^ — > z/ e ) and Viy^ — > v e ) — Vip^ — > v e ) 
since the neutrino beams at MiniBooNE, JHF-Kamioka, and MINOS, are 
muon or anti-muon neutrinos. We then calculate v e disappearance probabil- 
ity for Daya Bay baseline and energy, comparing our S-Matrix theory to the 
formula used by Daya Bay, Double Chooz and RENO|19 s 



2 Transition Probability V{v^ — > u e ) and v e 
Disappearance 

In this section we review z/ M to v e oscillation probability derived from stan- 
dard S-martix theory and then compare the probability of v e disappearance 
derived from this theory to that used by the Double Chooz experimental 
project [211 [22]. 



2.1 Viy^ — > u e ) Derived Using S-Matrix Theory 

In this subsection we review the derivation of the probability of a muon neu- 
trino to convert to an electron neutrino, Viv^ — > u e ), using the notation of 
Ref[31]. We then make an estimate of the transition probabilities for sample 
accelerator and reactor experiments. Although at the present time no exper- 
iments for CPV are possible, this can serve as a basis for future experiments. 
In the next section we give somewhat more accurate calculations for CPV 
for the same set of experimental facilities. 

As in Refs[2"9~| I5T] we use the time evolution matrix, S(t,t Q ) to derive 
the transition probabilities. For neutrino oscillations the initial neutrino 
beam is emitted at time to, usually taken as to — 0, and the neutrino or 
converted neutrino is detected at baseline length = L at time=t. Since the 



neutrinos move with a velocity near that of the speed of light, at the end of 
our derivation we take t — to — > L, with the units c=l. 

Given the Hamiltonian, H(t), for neutrinos, the neutrino state at time = 
t is obtained from the state at time = to from the matrix, S(t,to), by 

\u(t)> = S(t,t )\u(t ) > (2) 

ij t S(t,t ) = H(t)S(t,t ). (3) 

Neutrinos (and antineutrinos) are produced as v a , where a is the flavor, 
a = e, fi, t. However, neutrinos of definite masses are u a , with a = 1,2,3. 
The two forms are connected by a 3 by 3 unitary transformation matrix, U: 
v a = Uu a , where z/ a , v a are 3x1 column vctors and U is given by (sinOij = s^-) 



C12C13 S12C13 si 3 e zScp 

l r ( -S12C23 - ci2S 2 3Si 3 e 4<5cp C12C23 - s 12 s 2 3Si 3 e lScp s 23 c 13 

Sl2S2 3 ~ Ci2C23Si3e* CP — C12S23 — Si2C23Si3e* CP C23C13 



similar to the CKM matrix for quarks. We use the best fit value [12] S23 = 
0.707. 0i3 and s 12 = 0.56, c 12 = 0.83 . We use s 13 = 0.19 and Si 3 =0.095, 
as discussed above, to determine the dependence of v^ — > v e conversion, and 
CPV on this parameter. We calculate the dependence of Viy^ — > u e ) and 
AV^f on 5cp, as discussed below. In the vacuum the S(t,to) is obtained 
from 

3 

S ab (t 7 t ) = J^^ 6 ^^'^ 11 ^ ■ ( 4 ) 

i=i 

Since neutrino beams in neutrino oscillation experiments travel through 
matter, and the main neutrino-matter is scattering from electrons, we must 
include potential, V = v2Gpn e , for neutrino electron scattering in the earth: 
where Gf is the universal weak interaction Fermi constant, and n e is the 
density of electrons in matter. Using the matter density p=3 gm/cc, the 
neutrino-matter potential is V = 1.13 x 10~ 13 eV. 

The transition probability Viy^ — ¥ v e ) is obtained from Si 2 : 

%-H/ e ) = l^^e^f + Jm^] 2 , (5) 

with 



Su - 


= c 23 fi - is 2 3ae- lScp A 


(6) 


a = 


= s 13 (A - s 2 2 cT) 


(7) 


5 -- 


= 5m 2 2 /(2£) 


(8) 


A = 


= <^ 2 3 /(2£) 


(9) 


A ~ /(*)/„* 


(10) 


i"a* = 


= / dt'a*(t')f(t') 
Jo 


(11) 


a(t) = 


= cosut — icos28sinu)t 


(12) 


/(*) = 


= e~ tAt 


(13) 


2lo -- 


= ^/8 2 + V 2 -25Vcos(26 l2 ) 


(14) 


P - 


= —isin26sinu>L 


(15) 


A = 


= A-(V + S)/2 


(16) 


in29 = 


5 
- S12C12— , 

UJ 


(17) 



where the neutrino mass differences are Sm 2 2 = 7.6 x 10 5 (eV) 2 and Sm 2 3 = 
2.4 x 10~ 3 (eU) 2 . Note that 5 < A, and £ -» L for t> „ ~ c. From Eqs.(5 to 
16): 

i?e[Si 2 ] = s 23 a[cos(AL + 5 C p)Im[I a *} - sin(AL + 5 CP )i2e[/ Q *] 
Im[Si2\ = —C23sin29sinuL — S23d[cos(AL + 5cp)Re[I a *] 

+sin(AL + 5 C p)Im[I a *]] (18) 

Using 5, uj <C A one can show that 

sinAL 

Re[I m \ ~ 



A 

T r T -, 1 — cosAL , , 

/m[/ a ,] ~ ^ . (19) 

From Eqs ffTglfTg]) we find 

V^^^Ue) ~ (c 2 3S 12 c 12 (Vw)smwL) 2 + 2(s23Si3) 2 (l-cosAL) (20) 

+2si 3 si 2 ci 2 s 23 c 23 ((5/u;)smc<;L 
(cos(AL + 5 C p)sinAL + sin(AL + S C p)(l - cosAL)) . 

We use Si3=.19 and .095 to show the effect of Si 3 . 
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Figure 1: The ordinate is Viv^ — > v e ) for MINOS(L=735 km), Mini- 

BooNE(L=500m), JHF-Kamioka(L=295 km), and CHOOZ(L=1.03 km). 
Solid curve for Si3=.19 and dashed curve for s 13 =.095. The curves are almost 
independent of Sqp- 



From Eqf l20p we obtain the results for Viv^ — > u e ) shown in Fig.l. These 
results can provide guidance for future experiments on CPV via z/ M <H- v e 
oscillation. Note that in Ref[8] Viv^ — > v e ) was calculated for the 295 km 
JHF-Kamioka project for E=0-2 GeV, and our calculation based on the the- 
ory developed in Refs.jSU |29], finds Viy^ — > v e ) is in agreement for E=. 4-1.0 
GeV with this earlier estimate. We calculated V{y^ — > u e ) for 5cp from 
-7r/2 to it/2, and the results are almost independent of Sep- The results for 
CHOOZ are shown in preparing for the following subsection on z/ e disappear- 
ance, even though Double Chooz, Daya Bay, and RENO projects have beams 
of v e rather than v^ neutrinos. 

3 z7 e Disappearance Derived Using S-Matrix 
Theory Compared to Daya Bay Evaluation 

In this subsection we derive v e disappearance, V{p e — > v e ), defined as 

V(u e -)• z/ e ) = 1 - V(u e -)• i/ M ) - V(y e ->• v T ) , (21) 

using the S-matrix method (see previous subsection), and compare it to the 
expression for V{p e — > v e ) used by the Daya Bay, Double Chooz, and RENO, 
which is (see, e.g., Refs [TTl l2Tj ) 

V VB {y e ^v e ) ~ l-4( Sl3 c 13 ) 2 S m 2 (^) (22) 

where A = Sm 2 3 /(2E), Eq(8), and Si3,ci3 = sin6\z,cos6\z. A third term 
with a factor of sin 2 (SL/ '2) [22] was dropped [TTj |2~TJ because 5m\ 2 <C 5m\ z 
and sin 2 (5L/2) <C sin 2 (AL/2) for Daya Bay baseline L=1.9 km and energy 
E=4 Mev. 

In the S-matrix method (see Ref|31j) the probability of v e oscillation to 
Up and v T is given by (see, e.g., Ref[31j) 

V SM {v e ^v,) = \S 21 \ 2 

V SM (u e ^u T ) = \S 31 \ 2 . (23) 

We take 5qp = 0, since | *S 12 1 2 is essentially independent of 5c p, so A = C 
(see Ref[29] for definition and proof). Therefore |5 , 2i| 2 =|5'i2(^ / — > — V)\ 2 , 
and |5 3 i| 2 =|5 , i2| 2 (V ->■ -V, c 23 -4 s 23 ,s 2 3 -4 -c 23 ). Using s| 3 = cf 3 = 1/2 
we find 

V SM {y e ^v e ) = 1 - [(AQ5sinooL/io) 2 + 2(si 3 ) 2 (l - cosAL)] , (24) 

with A = A + (V - 8)/2 and 2 Q = ^J5 2 + V 2 + 2<JVcos(20i 2 ) 

Fig. 2 for L=1.9km and Fig. 3 for L=10km are discussed below. 
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Figure 2: V{v e ->■ z/ e ). For s 13 = .15 (a) V VB and (b) P 5 ^ 

From Fig. 2 the ratios i21,i22, of 1 - V{p e ->■ z/ e ) for P DB (Eq(J221l) to 
P 5M (Eq(J23J) for sia = .15, and for s 13 = .15 for V DB and s 13 = .147 for 
V SM , for E ~ 4.0MeV and L=1.9 km are 



R\ 



m 



l-7^*( Sl3 = .15) 

" l-V SM (s 13 = .15) 

l-V SM (s 13 = .U7) 



1.04 
1.00, 



(25) 



which demonstrates that using the S-Matrix formulation for L=1.9 km and 
E~ 4.0 MeV one would extract s\ 3 = .147 from the data for which the older 
formalism finds S\ 3 = .15. This is a 2% correction. We use the notation E~ 
4.0 MeV as there is uncertainty in the antineutrino energy. 
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Figure 3: V(u e -> v e ). For s 13 = .15 (a) V VB and (b) V SM 



Fig. 3 is the same as Fig. 2, except the baseline is L=10 km, as future 
might use a longer baseline for a larger effect given si 3 . For E ~ 4.0MeV and 
L=10 km the ratios are 

«l- 1 'L°> = - 15) , = 1-47 

l-V 3M {s 13 = .15) 

*= l:p{Z::Z = L0 °- (26) 

Thus using the S-Matrix formulation for L=10 km and E~ 4.0 MeV one 
would extract s 13 = .095 from the data for which the older formalism finds 
s 13 = .15. This is a 35% correction. 

We have carried out similar calculations for the T2K project, with E=0.6 
GeV, L=295 km[18]. With both a larger L and larger E than Daya Bay, we 
find a correction of 2.4%. 

It is also important to note that our SM method gives V{p e — > v e ) ^ 1.0 
even for Si3=0. 



4 CP Violation AV^f 

In this section we shall extend the derivation of the transition probability 
'Piyp. -^ v e) °f the previous section to derive the CPV probability 

— \Sv2\ ~ I ^12 1 (27) 

with 512 defined in Eq(8) and 

S12 = c 23 /3-ts 23 ae l6cp A, (28) 

with ft = (3{V -¥ -V) and A = A(V ^ -V). For example (see Eqs(16,18)) 
2co = ^/5 2 + V 2 + 25Vcos{26 12 ) and A = A + (V - S)/2. Using conservation 
of probabiltiyjH], \A\ 2 = \A\ 2 . 

^V C J = 4(|/3| 2 - \P\ 2 ) - 2c 23 s 23 ai(Im[Pe iS ^A*] - Im[0 e <W])(29) 



From Eq fl29|) . the definitions in the previous section, defining s = sin(ooL), 
c = cos(uL) one finds 

2 —1 

APjf = cl 3 sj 2 cl 2 5 2 (^-^) + 2c 23 s 23 s 12 c 12 s 13 5(A-5si 2 ) (30) 

. . , s -r A — ojcosIQ s = T .A — ucos28 
(stnUcp[ — [c — cosAL) — =r- 1 (c — cosAL) — = )) 

v K u y ' A 2 -u 2 u A 2 -u 2 

s sinALiA — cocos29) + sinuLiuj + Acos29) 
-cos9cp{- ^ — - 2 (31) 



'LU A 2 - LU 2 

s sinAL(A — ujcos26) + sinuL{u + Acos29) s 

U! A 2 - U) 2 



(32) 



The results for AV^ e p for s 13 =.19 are shown in Fig. 4. Note that AV^f de- 
pends strongly on 5c p, which could lead to a measurement of this parameter. 
The large value of AP^ P for CHOOZ is promising for future experiments. 
AV^f is so small (from about 10~ 10 to 10~ 18 )for MiniBooNE, we do not 
show the results. Similar results for AV^f for si 3 =.095 are shown in Fig. 5. 
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Figure 4: The ordinate is AP(z/ M ->■ v e ) for MINOS(L=735 km), JHF- 
Kamioka(L=295 km), and CH00Z(L=1 km). sl3=.19, and a, b, c, d, e 

for 5 cp =tt/2, tt/4, 0.0, -vr/4, -tt/2 
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5 Conclusions 

We have estimated CP violation for a variety of experimental neutrino beam 
facilities, for values of the parameter S13 =0.19 and .095, and for Sep from 90 
to -90 degrees, since its value is not known. As our results show, the proba- 
bility Viy^ —7- u e ) is strongly dependent on S13 and is essentially independent 
of Sep (see Fig. 1), and therefore the measurement of Viy^ — \ u e ) should 
determine the value of the S13 parameter, as has been known for many years. 

Our new results for v e disappearance, as is being measured the Daya Bay, 
Double Chooz and RENO projects, however, make use of a different theoreti- 
cal formulation than that used by these projects to extract S13 from the data. 
We have shown that the recent result from the Daya Bay collaboration |24] 
with E~4 MeV and L=1.9 km, from which it was estimated that S13 ~ .15, 
by our analysis is sy& ~ .147, a 2% correction. This is small, but the goal of 
these projects is 1% accuracy for S13. For a baseline of L=10 km, with E~ 4 
MeV, we find S13 ~ .097 using the S-Matrix method, rather than .15, a 35% 
correction. Also, our SM method gives V{p e — >■ z/ e ) 7^ 1.0 even for Si 3 =0. 

The CP violation probability (CPV), A"P£, P , is strongly dependent on 
both of these important parameters. After the Double Chooz/Daya Bay/RENO 
determination of S13, both the JHF-Kamioka and Double Chooz projects 
might be able to determine the value of Sep, since for most of the values of 
Sep these projects would have nearly a 1% CPV, as shown in Figs. 4 and 
5. No experiments are possible now to test CPV via neutrino oscillations, 
since beams of both neutrino and antineutrino with the same flavor would be 
needed. However, in the future such beams might be available. Our results 
should help in planning such future experiments. 
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